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Abstract 

We explore the role of Majoron (J) emission in the supernova cooling process, as a 
source of the upper bound on neutrino-Majoron coupling. We show that the strongest 
upper bound on the coupling to Ve comes from the J process in the core of 

a supernova. We also find bounds on diagonal couplings of the Majoron to 
and on off-diagonal couplings in various regions of the parameter space. We 

discuss the evaluation of the cross-section for four-particle interactions {vn —> J J and 
nJ —> vJ). We show that these are typically dominated by three-particle sub-processes 
and do not give new independent constraints. 


1 Introduction 

The solar and atmospheric neutrino observations provide strong evidence in favor of neutrinos 
being massive. These experiments are sensitive only to mass squared differences [1], on the 
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other hand, the overall mass scale of neutrinos is strongly constrained by the Troitsk and 
Mainz experiments [2], Combining these pieces of information, we conclude that the masses 
of the three active neutrinos are very small. Among the plausible and economic models 
which are developed to give a tiny mass to neutrinos are Majoron models [3, 4], In these 
models, additional Higgs boson(s) are introduced such that their vacuum expectation values 
break the exact B—L symmetry of the model. The Goldstone boson associated with this 
symmetry breaking is called the Majoron particle, J. 

In principle, Majoron particles can interact with matter -electrons, nuclei and photons. 
However the cooling of red giant stars provides a strong bound on these interactions [5]. 
Hereafter, we will assume that Majorons can interact only with neutrinos. In the literature, 
two types of Majoron interaction have been studied: 


Ani = + glp^la^K) 

(1) 

Cint = hap^'^a^.{dJ)^p, 

(2) 


where J is the Majoron held, is a two-component representation of a neutrino of havor 
/5, Qap and hap are 3x3 coupling matrices. The matrix hap is Hermitian while Qap is a 
symmetric matrix. In the model [3], for a range of parameters, the interactions can be 
described by Eq. (1) (see the appendix of Ref. [6]). In this paper, we will use this form of 
the interaction however, as we will see later, in most cases our results apply for both forms. 
Also, we will not assume any special condition on the diagonal or off-diagonal elements of 
gap- Since we have chosen a general approach, our results apply to any massless scalar held 
that has an interaction of the form given by Eq. (1), independent of the underlying model 
for it. 

Majoron models are also interesting from the astrophysical point of view because, they 
provide the only mechanism for fast neutrino decay which has not yet been excluded ([7] 
and the references therein). The role of neutrino decay in the solar [8] and atmospheric [9] 
neutrino huxes has been extensively studied. The possibility of explaining the anomalies by 
pure neutrino decay is excluded. In Ref. [10], decay of solar neutrinos along with oscillation 
has been discussed and, for a normal hierarchical neutrino mass scheme, it has been found 
that 2 

192.1" = I - 3 X 10-= ■ 

In Ref. [11], the different aspects and consequences of decay of neutrinos emitted by 
supernovae, have been studied. Future supernova observations can provide strong bounds 
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on (or evidence for) neutrino decay and consequently on the Majoron coupling, provided 
that the uncertainties in supernova models are resolved. 

If Majorons are coupled to neutrinos strongly enough, they can show up in /5/5-decay 
experiments, changing the spectrum of the dual electrons. Non-observation of such an effect 
imposes a strong bound on the coupling constant [12]: 

|c/ee| < 3 X 10“^ 

Also, no sign of Majoron particles has been observed in the pion and kaon decays and 
therefore [13]: 

E Il7ezp < 3 X 10-^ and E IvP < 2-4 x 10"^ 

l=e,fi^T l=e^ii^T 

The strongest bounds on neutrino-Majoron coupling are obtained by studying the role 
of these particles in a supernova explosion. In fact, three types of bounds are obtained: 

i) Deleptonization: if the coupling, \gee\ is too large, Majoron emission can reduce the 
lepton number of the core of supernova via z/g —preventing the emission of an intense 
observable photon flux. In [14, 15, 16, 17] this effect has been studied; the result is 

\g,,\ < 2 X 10-^ 

This bound strongly depends on the details of the supernova explosion model. 

ii) Spectrum distortion: the production and absorption of the Majoron particle can affect 
the spectrum of the observed neutrino flux from a supernova explosion. This effect has been 
studied in Refs. [17, 15] and the result is 

\an\ = \Y.KiUina,^0\<io-‘- 

a,p 

This result suffers from the low statistics of the SN1987a data and can be improved by future 
supernova observations. 

iii) Energy loss: according to [18], the binding energy of a supernova core is Ei, = (1.5 — 
4.5) X 10^^ erg, which coincides very well with the energy emitted by SN1987a in 1-10 sec 
in the form of neutrinos. Hence the power carried away by any exotic particle such as 
Majoron cannot be larger than ~ 10®^ erg/sec. This imposes strong bounds on the coupling 
of Majorons. The effect of energy transfer due to Majoron emission has been studied in a 
number of papers [14, 15, 6, 19, 20, 21]. 

In the presence of matter effects, a number of three-point processes that are kinematically 
forbidden in vacuum become allowed. For example, neutrino decay becomes possible even 
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in the absence of neutrino masses. Also, neutrino annihilation into a massless Majoron, 
vv J, becomes kinematically allowed. The latter process has not been taken into account 
in previous studies. We will see that this is actually the dominant process contributing to 
energy loss in a supernova explosion. Previous studies must be reconsidered to take this 
effect into account. 

In addition, the previous papers either considered only g^e or studied the Majoron cou¬ 
plings collectively without attention to the interplay of diagonal and off-diagonal couplings. 
In this paper, we study the effect of Majoron emission in the cooling process of supernova 
core considering all the relevant processes. We hnd that even for very small values of cou¬ 
pling, interplay of different processes may change the neutrino densities inside the supernova, 
evading the bounds that would be valid without this effect. 

If the couplings are larger than some “lower” bounds, Majorons will be so strongly 
trapped inside the supernova that they cannot give rise to signihcant luminosity. Note that 
these “lower” bounds should be much larger than the limits at which Majorons start to 
become trapped. For such large values of coupling, Majoron production can completely 
change the density prohle of the core by transferring energy between different layers and by 
changing lepton numbers. In this paper we discuss Majoron decay and all other processes 
that prevent energy transfer by Majoron particles and derive the limits on coupling constants 
above which the produced Majoron cannot leave the core without undergoing scattering 
or decay. We do not attempt to calculate any “lower” bound on the coupling constants, 
because for large values of couplings, the density distributions inside the core need to be 
recalculated. However we evaluate four-point processes which become important for large 
values of coupling constants. In summary, there is an “upper” bound on coupling below which 
the rate of Majoron production is so low that it cannot affect the evolution of supernovae. 
The values of coupling above the “upper” bound up to a “lower” bound are not allowed. 
However, the values of coupling above the “lower” bound (which are also higher than the 
“upper” bound) are not forbidden by supernova cooling considerations because for such 
values of coupling, Majorons cannot escape the core freely. The forbidden range is then 
between the “upper” and “lower” bounds. In figure (1), we illustrate all the bounds on \gee\ 
to clarify the meaning of the “upper” and “lower” bounds. The shaded area is excluded by 
the supernova cooling process. 

This paper is organized as follows. In section 2, we calculate the cross-section of the 
relevant processes. In section 3, we briefly review the characteristics of the core. In section 
4, we derive the bounds on the coupling constants and the values above which the produced 
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Majoron will scatter before leaving the core. Conclusions are presented in section 5. 


2 Majoron interactions 

In this section we hrst introduce the Lagrangian. Then, in subsection 2.1, we derive the 
formulae for the neutrino propagator and the dispersion relation in the presence of matter. 
The interaction rates for different processes involving Majoron are derived in subsection 2.2. 

In the presence of matter, the Lagrangian of neutrinos can be written in the two- 
component formalism as 

c = - V„„)% - (3) 

where C = ia2, a and f3 are flavor indices, a = (1, —a) and map is the symmetric Ma- 
jorana mass matrix. The term ^‘^^Vap^p represents the matter effect. This term has a 
preferred frame, the frame of the supernova. In the flavor basis, V is a diagonal matrix; 
V = diag{Ve, Vjj, K) with 


Ve = V^ + Vc, V, = Vr = Idv, (4) 

where 

Cc = V2GpnB{Y, + nj, Vn = (5) 

Yi = {rii — ni)/nB and ub is the baryon density [22]. f In Eq. (5), the Ij^^-dependent terms 
are the result of neutrino-neutrino scattering. Since in the medium of interest (supernova 
core) and riu^ = np^,^ the corresponding Y parameters vanish and have been 

omitted from Eq. (5). In Ref. [24] it is shown that due to loop effects the values of and 
Vr are slightly different, however, the difference is negligible: Vjj — Er ~ 5 x lO-^Re [25]. In 
a typical supernova core, and 14 are of the order of 10 eV and 1 eV, respectively. 

Tt is shown in Ref. [23] that, if the neutrinos present in a medium are coherent superpositions of different 

flavor states, the off-diagonal elements of Vap can be nonzero. However inside the inner core the densities of 

Ve and Vfi are different and the densities of and Vt are very low and equal to the densities of and Dr, 

so the off-diagonal terms vanish. 

Tn section 3, we will see that these equalities are only approximately true [30]. 
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For the interaction term, we invoke the form of Eq. (1). But we note that the derivative 
form of the interaction in Eq. (2) can be rewritten using the equations of motion as 

Thus, for processes in which all of the involved states are on-shell (in particular, neutrino 
and Majoron decay, vv ^ J and vJ v) the two forms of interactions give the same results 
with the replacement 

9a(3 ^ TTlci'yh^j^'). ( 6 ) 

As we will see, the most important processes involve only on-shell particles. Therefore all of 
the bounds in this paper, apply for both derivative and pseudo-scalar forms of the interaction. 

Majoron is a Goldstone boson associated with the exact B—L symmetry so in vacuum it 
is massless. Inside the supernova core Majoron obtains a tiny effective mass, rUe//, due to 
elastic scattering off the background neutrinos. It can be shown that ~ \g\'^Ny/q where 
q is the typical momentum of the particles involved. For the values of coupling constants of 
order of the upper bounds the effective mass of Majoron is negligible {mljj/q -C K)- The 
effective mass can be considerable only if | 5 f| ~ 5 x 10“"^. 


2.1 The propagators and the dispersion relation 


After straightforward calculations, we hud 

^[(d ■ p - (7) 

cr,7 p • U ^ V-y 

(K(p)tl,(-P)) = '£c "''°2 ,, {%(.pH(-p)} (8) 

^ P ' <7 V a 

and 

($„(p)$J(-p)) = ^ (-p)) I G, (9) 

7 y' ^ \ * (3 

where all of the subscripts a, (3, 7 and a denote {e,/i,r}. As we will see, the diagrams in 
which these propagators are involved are important mainly when |p| ~ 17 so the effect of 17 
must be treated non-perturbatively. If the mass scale of neutrinos is high (mj, 3> vAm3), 
the masses are quasi-degenerate; niafi — In this case the formulae are simpler: 




iSaq 


ml - {p ■ a + Va){p ■ a - 17) 


(p-a + Va), 


( 10 ) 


proposed Tritium decay experiment, KATRIN [26], may be able to determine the mass scale. 
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(r (»)$!,( p)) = c 

(11) 

and 

{4„(p)4j(-p)> = -”"1/u - v \^'- 

mf,- {p-a + Vc,){p-a-Vc,) 

(12) 

Now let 

us hnd the dispersion relation. The Lagrangian (3) yields 



®a(p)^(-P ■ d - 14) = '^map^'^{p)C 

P 

(13) 

and 

(p- d - 14)<h«(p) = -^mapC^l- 

p 

(14) 

Expanding the states as 


®a(p) = 

■■ ^ Uaih,p)aaih,p) + vl{h,p)al{h,p) for which a-pua{h,p) 

h=-l,l 

^ k\p\Ua{h,p) 

we hnd 


(15) 


“■>('».?) = p„ + _ yj:m.,CvPh,p) 

(16) 

and 

(p"+ftipi K.)«„(A.p)=E (p„ _ 44; i2)«2(A.P). 

(17) 


To find the dispersion relation and energy eigenstates we should solve Eq. (17). Note that 
the dispersion relation depends on helicity. 

For m?/p <^V^p, one can easily show that 


~ p - ^ ^ ( 18 ) 

and that the mixing among the flavors is of the order of m‘^/ 2 p{Vi 3 — Va) 1 which can be 
neglected. 

2.2 The relevant decays and interactions 

In this subsection we hrst discuss the processes that produce Majorons, then we study those 
that annihilate or scatter them. For illustrative reasons, in the following discussions, we 
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ignore mixing {i.e., off-diagonal terms in both coupling and mass matrix) and we denote 
coupling, mass and effective potential by g, m and V, neglecting their flavor indices. In the 
cases that generalization is not straightforward, we will discuss the relevant steps. Before 
beginning the detailed analysis, we should discuss an important conceptual point. As we 
see in Eq. (18) the dispersion relation for neutrinos inside supernova is different from that 
in vacuum and hence some reactions that are kinematically forbidden in vacuum, can take 
place in the supernova core. As we will see the decay i? ^ iz-h J and the interaction z/z/ —> J 
(or V ^ VJ and vv ^ J depending on the sign of V) are kinematically allowed. 

In addition to these three-point interactions, there are other interactions that produce 
Majorons: 

• u + u ^ J + J and p + p ^ J + J] 

• v + v ^ J + J. 

As we will see the effect of the four-point interactions is negligible. 


2.2.1 u ^ u + J or u ^ u + J 


In medium, if V is negative (positive), the decay u ^ u + J (z/ —z/ -|- J) is possible. Let us 
suppose E < 0, then, without loss of generality, we can write 


Pp = {Pi - V,^,^,Pi) Pu = {pf + V,p/sin6',0,p/cos6') 


where we have neglected corrections of order of m^/pi V. Energy-momentum conservation 
implies that, 

Pj = {Pi -Pf - 2E, -pf sinO, 0,pi -pfCos9). 


Recalling that we have neglected the effective mass of Majoron, the process P ^ u + J 
is kinematically allowed if and only if pj = 0 and all the zeroth components of the four- 
momenta are positive. Pj = 0 implies 


1 — cos 9 = 


AV‘^-2V{pi-pf) 

‘^PiPf 


(19) 


For \V\ < Pf < Pi, the above equation can be satished with all of the energies positive. This 
means that the process is kinematically allowed. 

Restoring flavor indices, it can be shown that ior Vp + Va < Q the rate oi Pa ^ I'p + J is 
given by 

dr _ \gag?Pi-Pf,,, , ^ 


( 20 ) 



where pi and pf are the momenta of the initial and hnal neutrinos and pf extends from 
Max(||f4 + V^l, —Vis) to Pi- In the equation, we have also included the Fermi factor 

FpPf) = (1 - 

e T +1 

which reflects the fact that inside the supernova some states have already been occupied by 
neutrinos. 

Similarly, for 14 +> 0, the process z/q, —z //3 +J can take place. The decay rate is given 
by Eq. (20) replacing (pf) with {pf). The range of pf extends from Max(||I4 + I4|, Vp) 
to Pi. 


2.2.2 u + u ^ J or z/+z/— 

In vacuum, the processes u + or 9 + 9^J are not kinematically allowed. However 

in medium, where V is negative (positive) the process u + u ^ J { 9 + 9 ^ J) can occur. 
Let us suppose H < 0 and study the possibility of z/ + z/ —J. Without loss of generality, we 
can write the four-momenta of the initial neutrinos as 


Pi = {Pi + 0 , 0 ,pi) p2 = {p2 + V,p2sm9,0,p2cos9), 

for which pi + V and P 2 + V are both positive. Energy-momentum conservation implies 


Pj = {Pi +P2 + 2V, p2 sin 9,0 ,Pi+P 2 cos 9 ). 


Recalling that we have neglected the effective mass of Majoron, the process z/(pi) + z/(p 2 ) 
J{Pj) will be kinematically allowed if and only if pj = 0 or 


1 — cos 9 


2V^ + 2V(pi+P2) 
P1P2 


( 21 ) 


which can be satished for pi,p 2 > |V|. 

Neglecting V^/p^ effects, for + Vp < 0, it can be shown that the cross section of 
^a^l 3 ^ J is given by 


a = 


4 :pIpI\vi - V 2 


(Pi +P 2 )|V„ + Vp\S{cos9 


cos 6 * 0 ) 


( 22 ) 


where pi and p 2 are the momenta of the two initial particles, 9 is the angle between them 
and cos 6*0 = 1 — _ 

piP2 

Similarly, it can be shown that for Va + Vp > f), instead of z/^ + z//? —> J, the process 
t'a + Z //3 —J can take place with the cross section again given by Eq. (22). 
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2.2.3 The process z/ + z/ —> J + J 

For reasons that will become clear in a moment, we analyze v and v as wave packets rather 
than as plane waves. Let us ignore the neutrino mass for simplicity. Then, calculating 
diagram (c) in hgure 2, we hnd 

X ///(P2)9yft)92i^x (23) 

- (9 + r)=]) 

xa2i^{pi)f{pi)d^Pid^P2 + {ki ^ /ca), 

where ki and k 2 are the momenta of the Majorons and / f{p 2 )\P 2 )d^P 2 and //(pi)|pi)d^Pi 
represent the states of the neutrino and anti-neutrino, respectively. In Eq. (23), q = k 2 —p 2 
and we have considered the matter effects in the propagator: 

F^(gO) = 1 - l/[exp{{q^ - fx)/T) + 1] 


is the Fermi factor. 

For both positive and negative V, in the vicinity of {ki = pi, k 2 = P 2 ) and {ki = pi, k 2 = 
pi), there are poles which are non-integrable singularities. Without the wave packets, the 
total cross section would be divergent. Setting mu non-zero just shifts the pole a little bit 
and does not solve this problem. This is due to the fact that for negative (positive) V, the 
processes V ^ v + J and u + u ^ J { u ^ 9 + J and 9 + 9 ^ J) can take place on shell, 
so the singularity is indeed a physical one. Essentially, for E < 0 the reaction v9 ^ JJ can 
proceed in two steps, hrst, 9 ^ vJ and later, at a completely distant place z/z/ — J. In other 
words, the total cross section has two parts: i) a “connected” part; ii) a “disconnected” part 
which can be considered as two successive three-point processes. 

Let us now consider in more detail the relation between Eq. (23) and its component 
three-point processes. For dehniteness, we consider the case E < 0. We have explained that 
the reaction 9v ^ JJ contains a subprocess that factorizes as 
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where F^(g°) = 1 — l/(ea;p((g° — fi)/T) + 1) represents the matter effects. In Eq. (24), 
r represents the boundaries on time integrations and therefore it must be very large (he., 
r ~ 5/E). We have written the time boundaries explicitly to emphasize the causality 
conditions. Transferring the amplitude for Pz/ —JJ (Eq. (23)) from momenta pi, p 2 to 
coordinates Xi, 0 : 2 , it can be shown that for the region \x 2 — Xi\ > 2r ~ 10/E, these two 
correspond. Therefore, if the initial neutrino and anti-neutrino are localized at distance 
R > 10/E, their interaction rate can be calculated by Eq. (24) instead of Eq. (23). 

Consider v and P which are localized at distance R > 10/E far from each other. We 
have shown that their interaction cross section is given by So, this 

interaction can be considered as two subsequent processes. First P decays into Ji and z/. 
Then, the produced neutrino propagates a distance R and annihilates with the other v into 
J 2 . In other words, to calculate the interaction probability of two such states, we can consider 
P —z/J as an additional source for v and consequently the process vv J. This can be 
compared to the more familiar sources of neutrinos like electron capture, {Jluu){un\e~p'^). 

Of course the set of states that are localized at distance R > 10/E ~ Rcore far 
from one another is not a complete set. We should also consider the states which are closer 
and/or have overlap with each other. If we rewrite Eq. (23) in the x-coordinates, as we have 
done in Eq. (24), calculation of the amplitude of two states localized next to each other at 
distance R will be easier. For such two states, the integral for \xi — X 2 \ > 10/E vanishes 
(because of the specihc form of f{pi) and /(P 2 )), so we can restrict the integration over 
\xi — T 2 I to the interval (0, 10/E). For \q^ — V —\q\ \ 3> |E|/10, the amplitude for two states 
localized at i? < 10/E far from each other is equal to the amplitude for states with definite 
momenta, but for \q^ — V — \^ \ < |E|/10, the amplitude for the two localized states is much 
smaller. This is because in calculation of amplitude for two states with dehnite momenta, 
we encounter an integration which diverges for q^ — V — \^ —> 0 

but for two states which are localized next to each other the corresponding integration 
is which is hnite. The total cross section for neutrinos and anti¬ 

neutrinos localized next to each other is then given by an angular integral over the square of 
(23) in which the integration is over all angles except those for which — E — |^ | < |E|/10. 
Consider the special case that the sum of the momenta of v and P is zero. Setting the cutoff 
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equal to A|l/|/10, (A is an arbitrary number of order of one) for such two particles we obtain 


87^PlP2|^'l - V2\ 


PIP 2 

(Al^/10)2 


10 14 

+ T“T 


Since we have a preferred frame (the frame of the supernova), the total cross section is 
not Lorentz invariant and for each 1 / and P the total cross section must be calculated, 
independently. Now consider a pair of a neutrino and an anti-neutrino that make a general 
angle. Then Eqs. (19,21) show that in the vicinity of singularity the momentum flowing in 
the propagator is of the order of |E| (he., = |p 2 — ^ 2 ! ~ |h^|)- Therefore for scattering 

angles that {<f — V-q-a)/ [(g° — 1/)^ — ~ 1/E, the phase factor (/ d?kid?k25^{pi+P2 — 

^2 — ki)) is of the order of \V\‘^/p^. Thus the total cross section has no strong dependence 
on |E| and for general initial momenta the cross section can be estimated by Eq. (25). 

Here, for simplicity we have dropped the flavor indices but for the more general case the 
discussion is similar. 


2.2.4 The process u + u ^ J + J and u + u ^ J + J 

The discussion of uu ^ JJ and P + 9 ^ J -|- J can be carried out in a similar way. For 

quasi-degenerate neutrino masses, the amplitude for z/z/ —JJ (see diagram (d) in Fig. (2)) 
is given by, 

-(26) 

{2TiY^Ak\kl^ J J 

irn{m? + V^ + — ql + 2q- aV^) 

{w? - ql + (Ey - qy){w? - ql + (E^ -h qY) ^ 

f2{p2)u{p2)<Ypi<Yp2 + {ki ^ k2) + A, 

where / fi{pi)\pi)d^Pi and / f 2 iP 2 )\P 2 )d^P 2 represent the initial neutrino states, ki and /c 2 
are the momenta of the emitted Majorons and q = k 2 —p 2 - The term A summarizes all of the 
Fermi effects on the propagator. The amplitude for values of q which {q + E^)^ — gg — ~ 

PitP 2 ^ is negligible, and the main contribution to the cross section comes from the 

small solid angle (~ E^/piP 2 ) for which (g -|- E)^ — gg — ~ E^. 

First, let us discuss the process JJ. In general, for 7 = /i, r, there are singularities 

which correspond to an on-shell z/^ or Note that, if pi and p 2 are parallel or make an 
angle smaller than ~ |Ee/E^|, the singularities disappear. As for the case z/z/ —JJ, we can 
discuss that if the initial states are localized at distance R > 10/V/ far from each other, 
the process z/z/ —JJ will be equivalent to two successive processes {i'^(T)J\ve) and then 
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(J|z/eZ/^(^)). This yields a cutoff of V^/IQ for calculating the 4-point total cross-section. Note 
that although z/g —^ J + r'/i is kinematically allowed {V^ < 14 ), T{ue Jv^) is suppressed by 
{rn/pu^Y practice, will not have any signihcant effect. 

For 7 = e, there is no singularity (except for the case that one of the Majorons is soft 
and z/gZ/g —J is kinematically possible) and therefore no cutoff is needed. The total cross 
section can be estimated as 


^ JJ) = 

Similarly, 


|Ul - U2|(27r)2piP2 ' 

^(t)) ^ t/T) 




|ui - U2|(27r)2piP2 

\9e^l9^l^l{T) + 9eT9T^i{T)\ )(y~/ 10 ^ ^ \9ee9e^l{T)\ ^ 


and 


^(t) ^ t/T) 


|ui - U2|(27r)2piP2 


+ 9%(t)? + ^"\9eti{T)f) 


(27) 

(28) 


(29) 


with b, b', 6 ", a, a' and a" of order 1. In Ref. [20], <Jtot{^e^e JJ) has been calculated, 
ignoring V and the off-diagonal elements of the coupling matrix. The result agrees with our 
estimation in the sense that the term proportional to \9ee\‘^ is not suppressed by m. 

The total cross section for (z?aZ/^ — JJ) is equal to JJ) replacing V with 


2.2.5 The processes z/-|-J—>z/orz/-|-J—>z/: 


These processes are the opposite of anti-neutrino and neutrino decay and, hence, the kine- 
matical conditions are similar. \i + Vp is negative (positive) the process z/^J —( 
z/q, J —z//j) can take place with cross section 


a = 


( 2 '?’') |„ |2 |ha + hal 


'Ifi'o/sl 


F(-){p + q)5(cos6 — cos 6 *o) 
P 9 


(30) 


4pg|ui - U 2 I 

where p and q are the momenta of the initial neutrino and Majoron, respectively. 6 is the 
angle between the two initial states and 

{p + q)K + Vp\ 


cos 6*0 = —1 -|- 


pq 


FlL){p + q) is the Fermi factor for the hnal state. 
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2.2.6 The Majoron decay, J—>z/ + z/or J—>z/ + z/ 

The decay J ^ uu {J ^ vv) is the opposite of the interaction vv ^ J {99 —> J) and 
therefore the kinematics are similar. 

For Va + Vjs < Q {Va + Vg > 0) the Majoron can decay into {9^ + 9/^) and the 

decay rate up to a (|V^|/pi)^ correction is given by 

<ir = [’ F^iPl)Fs (P. - Pl)dPl (31) 

where p* and pf are the momenta of the Majoron and either of the hnal neutrinos, respec¬ 
tively. and are the Fermi factors reflecting the fact that in the core of the supernova 
some states have been occupied by already present neutrinos. 


2.2.7 The processes u + J ^ u + J and u + J ^ u + J 

The amplitude for + J 9^ + 3 has two singularities in the t-channel due to z/^ ex¬ 
change. Using a similar discussion to that in section 2.2.4, it can be shown that these 
singularities may be considered as two successive three-point interactions {9e\Jiy^){iy^J\iye) 
and This yields a cutoff ~ |Vjj|/10 around the singularity to determine the 

four-point interaction. In the case of head-on collision where the initial particles are within a 
small solid angle ~ {V/pY ^ 47r around cosd = —1, there will be another singularity in the 
s-channel which can be considered as {9eJ\9^^){9y\veJ)■ We recall that any discussion about 
z/^ applies to Vr as well, because these states are completely equivalent for the supernova 
evolution. The total cross-section for 9^3 can be evaluated as 


(27r)2|r;i - n2|piP2 


a\g% + gU\ 


m 


vpw 


+ i>l9t 


(32) 


where a ^ 6 1 and is the Fermi-blocking factor for the final neutrino. A similar 

discussion holds for z/gT —^ z/gJ, and the corresponding cross-section is also of the form of 
Eq. (32). 

The processes h'^ + J^9e + J,Ue + J^9^ + J,9fj_ + J^h'e + J and 9f, + J ^ J 
also have singularities in the t-channel due to z/^-exchange and can be considered as two 
successive three-point processes. Following the same discussion as in sections 2.2.3 and 
2.2.4, we use the cutoff ~ Vji/lO to evaluate the cross section for the four-point interactions. 
The cross-sections of these processes have the form 


(27r)2|t;i - V2\piP2 


(^\gefj,g^^ T gerg^^l 
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KVio 




(33) 
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where a ~ 6 ~ 1 and is the Fermi-blocking factor for the hnal neutrinos. The processes 
> VeJ and can also have singularities in the s-channel only if the initial 

particles are almost parallel, i.e., if their relative angle resides within a small solid angle 
~ (V/pY 4:71 around 180°. We can safely neglect such states. 

For the process ^ J, there is no singularity and it is straightforward to show that 
the cross section is of the form, 

+ «?..! +''IviT (34) 


2.2.8 The processes z/ + J —> z/ + J and z/ + J —> z/ + J 


In general, the process z/ + J ^ z/ + J has a singularity in the ^-cha nn el. With a similar 
discussion as in section 2.2.3, we can show that this singularity can be evaluated as two 
successive three-point interactions {J\vv){yv\J) resulting in a cutoff of the order of V/IQ for 
evaluation of the four-point interactions. Using this cutoff, the cross section is of the order 
of 


w 


In 


(JP^\ pF 


(35) 


{277f\Vi -V2\piP2 V^VlOOy 
where is the Fermi-blocking factor for the hnal neutrino. 

If the initial particles undergo a head-on collision {i.e., they are within a small solid 
angle ~ {V/pY around 180°) there will be another singularity in the s-channel which can be 
considered as (z/J|p) (P|z/J). The process u + J ^ 9 + J has one singularity which can be 
evaluated as (P| Jz/)(Jz/|p). Again the cross section is of the form of Eq. (35). 


3 Supernova core without Majorons 

The dynamics of a supernova explosion is described in a number of articles and books {e.g., 
[22]). Here we only review the aspects of the supernova explosion which are relevant for our 
calculations. 

Very massive stars {A4 > S/Aq), at the end of their lifetime, develop a degenerate core 
with a mass around l.SAf© made up of iron-group elements. As the outer layer burns, 
it deposits more iron that adds to the mass of the core. Eventually the core reaches its 
Chandrasekhar limit, at which the Fermi-pressure of the electron gas inside the core cannot 
support the gravitational pressure, and the star collapses. The collapse forces nuclei to absorb 
the electrons via e~ + p ^ n + At the early stages, the produced z/g can escape from the 
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core but, eventually, the core becomes so dense that even neutrinos are trapped. The layer 
beyond which neutrinos can escape without scattering is called the “neutrino-sphere”. 

As the density of the central core reaches nuclear density (p ~ 3 x 10^^ g/cm^), a shock 
wave builds up which propagates outwards. We will refer to the pre-shock stage as the infall 
stage. This stage takes only around 0.1 sec. As the shock wave reaches the neutrino-sphere, 
it dissociates the heavy nuclei. The dissociation has three different results: 

1. It consumes the energy of the shock, so that the shock eventually stalls; 

2. It allows neutrinos to escape more easily; 

3. It liberates protons that interact with the electrons present in the star {e~+p n+Ue), 
giving rise to the famous “prompt z/g burst”. The prompt z/g burst deleptonizes the 
star but carries only a few percent of the total energy. 

The stalled shock should regain its energy. Otherwise, it cannot propagate further and give 
rise to the spectacular fireworks. According to the models, this energy is provided by z/g 
diffusing from the inner core to outside. The density of z/g inside the inner core is very 
high. The corresponding Fermi energy is ~ 200 MeV while the temperature is only around 
10 MeV. At the beginning the temperature of the neutrino-sphere is around 20 MeV. So 
the diffused neutrinos leave their energy as they travel outside, warming up the core. This 
energy can revive the shock. (In fact, this mechanism is controversial [22], but we will not 
use the shock revival mechanism for our calculations. Most of our calculations are related 
to the inner core, which is free of these controversies.) The temperature in the outer core 
increases to 40 MeV; actually, the outer core and the neutrino-sphere become warmer than 
the center. At the onter core, nentrinos of each type (z/g, z/g, z/^, z/^, and 9^) are present. 
These neutrinos escape the star and deplete its binding energy {Ef, = (1.5 — 4.5) x 10®^ erg 
118]). 

Two kinds of “upper” bounds can be imposed on the nentrino-Majoron conplings by 
stndying snpernova evolntion: 

1) If the coupling constant is too large, the process z/g —J -|- z/g, during the infall stage, 
deleptonizes the core and according to the models [27] a successfnl explosion cannot occur. 
This bound has been correctly studied in [14, 15] and the resnlt is Qee ~ 2 X 10-6. 

2) If the conpling is non-zero, Majorons can be produced inside the inner core and can 
escape freely from the star, depleting the binding energy. The observed nentrino pulse from 
SN1987a coincides with that predicted by current snpernova models. This means that the 
energy carried away by Majorons (or any other exotic particles) should be smaller than 
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t(sec) 

W (eV) 

Vp = Vr (eV) 

0 

2.3 

-11.7 

0.5 

1 

-12.3 

1 

-0.3 

-12.8 

1.5 

-1 

-13.1 


Table 1: The values of the effective potentials at different instants after bounce without 
Majoron production. 


the binding energy. The Majoron luminosity, £j, as large as 10^^ erg/sec could signihcantly 
affect the neutrino pulse. Here, we will take £j < 3 x 10^^ erg/sec as a conservative maximum 
allowed value. This gives an upper bound on the coupling constants. 

If the coupling of Majorons is larger than a “lower bound”, the Majorons will be trapped 
so strongly that their luminosity will be small. We will discuss this case later. 

Let us review the characteristics of the core. The inner core {R < Rinner ~ 10 km) to 
a good approximation is homogeneous. The density in the inner core is around 5 x 10^^ 
g/cm^. The distributions of all types of neutrinos follow the Fermi-Dirac formula with 
Tinner ~ 10 — 30 MeV and different chemical potentials [28, 29]. As mentioned earlier, the 
chemical potential for is around 200 MeV. So, inside the inner core, is degenerate while 
the density of is negligible = —200 MeV). The suppression of the density of 

Ve is due to absorption on electrons. In the hrst approximation, the chemical potentials for 
and are equal to zero. In Ref. [30], it is shown that, because the interactions of 
and Vr with matter are slightly stronger than the interactions of and Vr-, their chemical 
potentials become nonzero: fiu^/T = Hur/T — bT/rrip < 1. We will neglect and in 
our analysis. In fact the large uncertainty in the determination of the temperature affects 
our results more dramatically. The presence of /i in supernova can break the equivalence of 
Vp and Vr- However, we neglect this effect and treat and Vr in exactly the same way. In 
Table 1, we show the values of and Vp (= W) at different instants after the bounce inside 
the inner core. The values of W and are taken from Ref. [28]. 

Outside the inner core, Rinner ~ 10 km < R < Rout ~ 15 km, the density of z/g 
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is much lower, ~ 1, but instead the density of is higher than in the inner core. 

In fact, in the outer core {Rinner < R < Rout), thermal equilibrium for neutrinos is only 
an approximation. To evaluate the role of the outer core in Majoron production, we set 

= fJ'Uy = = 0- The density in the outer core drops from 5 x 10^^ g/cm^ to 5 x 10^^ 

g/cm^. The temperature in the outer core drops abruptly [28] such that T{R = Rinner) = 35 
MeV while T{R = Rout) ~ 2 MeV. 

Different models predict different values for parameters; e.g., the predictions of different 
classes of models for Tinner vary from 10 MeV to 30 MeV [29]. Moreover the production of 
Majorons can distort the density distributions. Considering these uncertainties, the simpli- 
hed model that we have invoked is justihed. With this approach, we will be able to examine 
the prediction of all models for the Majoron luminosity. 


4 Bounds on coupling constants 

In this section we explore the role of Majorons in the cooling of the supernova core. In 
subsection 4.1, we derive an upper bound on \gee\ assuming the produced Majorons leave the 
core without being trapped. In subsection 4.2, we derive upper bounds on g^r, g^ie and 
gre, again assuming Majorons leave the core immediately after production. In subsection 
4.3, we show that for the couplings lower than the bounds we have derived, the four-point 
interactions are negligible. In subsection 4.4, we derive the limits above which Majorons 
become trapped. 


4.1 Bounds on l^fgel 

As represented in Table 1, immediately after the bounce, 14 is positive, but eventually 14 
decreases and becomes negative, while Vjj and 14 are negative from the beginning. As long 
as 14 > 0, the interactions z/g —> he + 4 and Ve + J are kinematically allowed but the 

latter is suppressed by a factor of (m/p)^ ~ 10“^®. So we will consider only the interaction 

z/g ^ he -f J. 

This interaction depletes the energy of the core at a rate 

X W^Rlner)- (36) 
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We should note that this interaction not only carries energy away but also deleptonizes the 


core. 


dYr 


= - 2 rw = -2 


gleVe 




(37) 


dt Stt 

where we have used the fact that We know that the core is in /^-equilibrium. Since 

the rate of the /5-interaction is faster than T (rate of /3-interactions/ T 

and at equilibrium the density of electrons is one order of magnitnde larger than that of 
neutrinos, we expect that the densities of the neutrinos are not affected by Majoron produc¬ 
tion. In other words, the Fermi energy, and Y/e ^ire still given by Ref. [28]. However, 
deleptonization by Majoron emission can affect V/ dramatically because different terms in 
Ve oc (31 l + Yiy — l)/2 cancel each other (Y/ Yl — 0.3). Therefore in the presence of 
Majoron emission, V/ vanishes faster. Let us evaluate the maximum energy that can be 
carried away by Majorons throngh z/e —zzg -f J in the stage when V/ is positive. To have an 
estimation, we can approximate 


dt 


= -bV, 


(38) 


where 


b = V2^GF^\g, 

Stt Mat 


W 

J- ly 


and a reflects the deleptonization effect without Majoron emission. According to Table 1, 
a ~ 2.6 eV/sec. If we neglect the variation of Y/, p and a with time, we conclnde that 


Veit) = (14(0) + 

so that, after ti = (1/&) x ln(14(0)6/a-|-l), 14 vanishes. The energy carried away by Majorons 
up to ti can be approximated as 


Eve<o — 


12(27r) 


X 4/37ri?; 


3 

inner 


^^^14(0) 14(0)6-ha, 

b 62 " a 


(39) 


For gee 10 Ev^<o converges to 4 x 10^^ erg. Increasing gee increases Cj, bnt on the other 
hand, 14 vanishes in a shorter period. It is easy to show that, for any valne of gee, 


Eve<o < 4 X 10^^ erg < E^ 


Therefore the energy loss at this stage does not affect star’s evolntion and hence we do not 
obtain any bonnd. 

As shown in Table 1, about one second after the core bounce 14 turns negative. As we 
discussed earlier, in the presence of neutrino decay 14 changes its sign even faster. In a 
medium with negative 14, the decay z/g —> zzg -|- J is not kinematically allowed and instead 
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i>e —z/e + J can take place. However, we know that, in the inner core, the density of electron 
antineutrinos is quite low ~ —200 MeV while T ~ 10 MeV) so this interaction will not 
have any role in the cooling of the inner core. In such a medium, energy will be carried away 
by process 

+ j- (40) 


In previous literature the possibility of this interaction was not discussed. The interaction 
(40), diminishes the lepton number by two units. Again we see that and will not 
be considerably affected by this process, but that Vg will decrease faster. In contrast to the 
previous case, a faster decrease of lA is a positive feedback for the process and leads to the 
energy depletion. The energy carried away from the inner core via the process in Eq. (40) 
is now 

X (41) 

To evaluate a conservative upper bound on \gee\i we set |I4| equal to 0.3 eV, jjiy^ = 200 
MeV and i? = 10 km then. 


Cj = 2\g^e? X 10 ^^® X 


Ri, 




:)(: 


P'l- 


erg 

sec 


10 km^ ^0.3 eV^^200 MeV^ 

Around one second after the core bounce, the total neutrino luminosity, Cy, is about 5 x 10^^ 
erg/sec. So, the condition Cj < 3 x 10^^ erg/sec yields the conservative bound. 


\gee\ < 4 X 10 \ 


R, 


inner \ — i 






1-2 


(42) 


10 km' 0.3 eV '200 MeV' 

In Ref. [14], a bound on \gee\ is obtained by studying the energy loss via Ve ^ + J 

which mainly takes place in the outer core, Rinner — 10 km < R < Rout — 20 km. 
The result is C{ve Ve + J) = few x 10®"^ Ifi'eeP erg/sec. So the conservative bound 
C{ve Ve + J) < 3 X 10^^ erg/sec implies \gee\ < 4 x 10“®. The bound in Eq. (42) is one 
order of magnitude stronger because the total number of in the inner core is very high. 
In Ref. [15], a bound is imposed due to the processes z/ + z/ J + J and v ^ u + J {y 
denotes both neutrino and antineutrino). However the energy carried away is overestimated 
due to an improper treatment of the three-point subprocesses. We will elaborate on the 
p + p ^ J + J process in section 4.2. 


4.2 Bounds on \g^a\ and l^r^l 

In this subsection we discuss the processes involving and/or These processes include 

(®) ^/4,T T ^ J-i ^/4,T T ^ J 
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and 


The process ^ J + ^e can take place only in the outer core where electron neutrinos are 
not degenerate. Both processes (a) and (b) can distort the distribution of matter inside the 
star. However, that calculation is beyond the scope of this paper. But we can argue that it 
is a good approximation, for the purpose of computing upper bounds, to use distributions 
with vanishing chemical potentials for and [28]. For simplicity we rotate (z/^, vP) to a 
basis such that = 0. Note that since the chemical potential is diagonal and = fd-, it 
will be invariant under this rotation. In the new basis, we can write, for the inner core 


dn,. 


dn,-. 


= 2 [Rate(z/^ + J) — R-ate(z/^z/^ —J)] — Rate(z/^z/e —> J), 


(43) 


dt dt 

where we have neglected vv ^ JJ interactions. The sum of the chemical potentials for 
and must be zero, /i = = —/ip,,, therefore 

dTT p^dp 


= 


An p^dp 


(44) 


(27r)3e^ + l- 

We expect that for small values of \gap\i the chemical potential remains small. Let us suppose 
|/i/T| 1 to solve the equation (43), then we can determine whether this assumption is 

valid or not. For |/i/T| 1, 

47rT3 


- 


(2ir) 


r(1.8 + 1.64>i/r), n„ ~ (4irr7(2ir)'’)(1.8 - 1.64;i/r) 


and we can rewrite the right hand side of Eq. (43) as 

-{0.12 - 3.28/i/T - (0.34 + 0.25/i/T) 




lde/i| 

'|V; + Re| 






}. 


(45) 


2 ( 2 ir)s 

where is the chemical potential of the electron-neutrinos. Inside the supernova core, 
neutrinos and matter are in thermal equilibrium and since the energy density of matter is 
much higher, we expect that the rate of thermal change due to these processes is small: 

dT E^^driy^/dt dn„ 


'Tdt' 




n,y^dt 


Efe/volume 

So, driy^/dt — dup^/dt ~ x 1.6Ad{p/T)/dt. On the other hand, for this estimation we 

can neglect the variation in V)j ~ ^/^GpUBiXe — l)/2. Also, since the density of z/g is much 
higher than that of z/^, we can neglect the variation of Therefore, Eq. (45) tells us that 
p/T converges to 


(0.12-0.34 


\9e)i\ 

"\Ve + V^\ 


\9^l^l 




)/(3.28 + 0.25 


1 de/ll 

^\Ve + V^\ 



? I^.l 
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Now, it is easy to show that for x (/i^^/T^) < 37, |/i/T| remains small regardless 

of the values of \geij\ or themselves. For \geij\ > Q\gtM^l\T/ Ih/^l diverges to values 
larger than 1 and the above analysis is no longer correct (remember that we had assumed 
|/i/T| 1). In this case, will disappear after ~ ( ^looi but on the other 

hand, the density of 9^ will increase (the chemical potential becomes negative) and this calls 
for recalculation of the density distributions. We can make a similar discussion for Ur. Let 
us suppose \ge^\ < and {gerl < ^IgrAT /and continue from here. 

Now let us evaluate the Majoron luminosity using the distributions given in [28]. Ne¬ 
glecting the Majoron emission from the outer core we hnd, 

2 

+ ^-^TA^neAgcAWVo. + vj,|) (46) 

and 

1 3 R^ 

C{9^ ^ upj) ^ + f^/3l), (47) 

where by a and f3 we denote g, or r. On the other hand, 

d3 

C{iy^ + z/e ^ J) ~ + W|)(0.2/i3 T„). (48) 


Note that even if g^e is large, the process ( z/q, —z/g + J), in the inner core, is suppressed by 
a factor of exp{{Tinner — RuJ/Tinner) bccausc inside the star, z/g is degenerate. 

Then, the requirement £ < 3 x 10^^ erg/sec implies 


/ ^ ~ 8 X 10-b 

a,g£fi,T 


,10 eV.1,20 MeV, 
\2 -' 
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(49) 


and 


n --;-rr _ 7 /J-U Km, 3 , 

+ IdreA < 5 X 10 - )2( 

■Ttir. 


10 km, 3,200 MeV,3,20 MeV, i TO eV,i 


2(- )-2[— -(50) 

^ inner ^/i 

We emphasize again that the above results are valid only assuming that \gei_igue/gt^iiTA < 37 
and \geT-gue/grrTA < 37. Otherwise the annihilation will stall because is depleted. 
Meanwhile, the energy carried away due to z/^(t-)- annihilation is of the order of 


£(^e 4“ ^ "7) ( 


1 ^. 


IOOtt 


10"^® erg Eb ~ 10^'* erg 


53 


So, £(z/g -f ^ J) does not impose any bound on \gefj,g)\- On the other hand, since the 

density of 9^g) grows, the process Cip^g-) —z/^(r)<7) will even become intensihed and we 
expect that still Eq. (49) will be a conservative bound. For {\geA‘^l\9^^AA ^ {rI^/TA > 37, 


22 



the upper bound on \getj\ is imposed by z/^-decay in the outer core. Using the distributions 
in Ref. [28], we can show that 

J + Ve) = few X X 10®^ crg/scc (51) 

which implies 

\geA, Ifl'erl < feW X 10“®. (52) 

In Figs. ( 1 ) and ( 2 ), all these bounds are schematically depicted for Tinner = 10 MeV and 
Tinner = 20 MeV, respectively. The shadowed area represents the range of parameters for 
which Cj < 3x 10 ^^ erg/sec. As shown in Fig. (4), for T = 20 MeV, the process + ^ J 

does not impose any bound on because, for any value of \geij,\ smaller than \^\g^^\T / 
(where is below its upper bound) it cannot give rise to a Majoron luminosity larger 
than the allowed value. 


4.3 Four-point interactions 


In this subsection we discuss the processes u + u ^ J + J and P + h —> J + J. As discussed 
in sections 2.2.3 and 2.2.4, we consider only the intrinsically connected contributions, the 
effects of three-particle snb-processes subtracted. Using the distributions in Ref. [28] and 
the formulae we have found in snbsection 2.2.3, we obtain 
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r^Ve inner 
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and 


T® dvr 

In the above equations, a runs over {e,/i,r}. Using the distributions in Ref. [28] and the 
formulae we found in section 2.2.4, we obtain 

1 ( TTh^ ^ 47r 

C{Ue + Ve ^ J + J) = y 1 / 2 + ^Ifl'eel J ( 

C{l'e + J + J) = (56) 


and 


+ 9er9r,ir)\^^ + b'\9ee9e,ir)\ 

U(z/q, 1^/3 —> J J) ^ + p/3 —J T t/) = 


(57) 
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3 (27r)6 

where m is the neutrino mass for quasi-degenerate mass schemes. If \gai 3 \ is smaller than the 
“upper” bounds in Eqs. (42, 49, 52) the above luminosities are negligible. These luminosities 
become non-negligible only if the couplings are larger than 10“^ so they do not change the 
“upper” bounds. Eqs. (53-56) depend on combinations of couplings so in general it is 
rather difficult to compare them with the Majoron luminosity due to three-point processes 
(h e., C{i?a J^p)) J))- Let us suppose that all elements are zero except for 

a particular g^p. Then, Eqs. (46,47,54,56) imply that for < 5 x 10“^ the three-point 
processes are dominant. If all couplings, but \gee\ are zero Eqs. (41) and (55) show that as 
long as \gee\ < 10“^, z/eZ^e —J is dominant. Comparing Eqs. (53) and (55) with Eq. (51) 
reveals that as long as \ge^l\ < 10“"^ the process z/^ —z/gJ dominates. We note that, for 
coupling constants of the order of the “lower” bound (for which the produced Majorons are 
trapped), the four-point processes can play a signihcant role. 


4.4 Majoron decay and scattering 

So far we have assumed that Majorons leave the star without undergoing any interaction or 
decay. Now we discuss the validity of this assumption. First, let us discuss the possibility 
of decay. (Note that although Majorons are massless particles, in a medium such as a 
supernova, in principle, they can decay.) For a,P E 

r(J(g) ^ z/, + up) = - 0.27), (58) 

87r 

where 0.8 and 0.27 correspond to q/T = 10 and q/T = 0.1, respectively. So, the Majorons 
decay before leaving the core (T > 1/77), only if 

\gap\ ~ lo-l (59) 

Because of degeneracy of the inner core, only the energetic Majorons {\Ej — fiue\/T ~ 1) 
can decay into electron neutrinos (see Eq. (31)). It can be shown that 

. o ^ ,1 T\gee\‘^\Ve\ 

T[J{q > 2/ij.J ^ z/e + z/e]--- 

and 

^ . \9ee?\Ve + V,\ 

r[J{q > ^ z/e -7 z/„]---. 
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If 


> 7 X 10 ® and/or \gee\ > 5 x 10 


(60) 


Majorons that are produced in the center will decay before leaving the core. Note that, even 
beyond the neutrinosphere as long as \V\ is large enough, Majoron decay can take place. 
Now let us examine the interaction effect. For low values of coupling constants, the dominant 
interactions are {u + J ^ 9) with the mean free path 

+ J^9,) = (61) 

47r q 




M in + vM 

Stt q 


+ 1 


(62) 


I {l'i3 + J —* l^a) — 




Stt 


-m+v^\) 


Te''/^(g/T + ln2 - ^(6*?/^ + 1) 


.q/T 


(63) 


and 


I (^/3 + "7 ^ «^e) — —X-- —{\Vf3+Ve\), 


(64) 
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where q is the energy of J and a and (3 are either fi or r. The requirement > R~^ implies 
that 
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and 


V2|^rMl ~ 4 X 10 

In the last case, the bound is derived for g = 10 MeV, T = 10 MeV and |V/| = 10 eV. 
Note that these bounds are derived for the parameters inside the inner core. In the outer 
core is much smaller and therefore the mean free path is much larger, i.e., in the outer 
core neutrinos can escape more easily. Apparently if the coupling constants are smaller 
than the bounds in Eqs. (42,49,50), Majorons will leave the star core before undergoing any 
interaction. For 6|g^^|(T//i„J < \ge^\ ~ 2 x 10“®, although the Majorons produced in the 
outer core escape immediately (recalling that the bound in Eq. (52) is extracted by studying 
the process —> z/g + J in the outer core), the interaction of Majoron particles with z/g in 
the inner core is not negligible. 

Eor larger values of coupling, Majorons may become trapped or decay before leaving the star 
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and the energy transfer by Majoron emission will become harder, but this does not mean 
that Majoron production does not affect the supernova evolution. To calculate the exact 
effect and to extract lower bounds on coupling constants, one needs to revisit the matter 
distribution and its time evolution including the effect of energy transfer by Majorons. That 
is beyond the scope of this paper. Here, we have discussed only the dominant interaction 
modes for larger values of the coupling constants. We recall that for \g\ ~ 5 X 10“^ the 
effective mass of the Majoron becomes non-negligible. 


5 Conclusions and discussions 

We have explored the energy loss from the inner core due to emission of Majorons. We have 
found that at early instants after the shock bounce (t ~ 1 sec) when W is positive, although 
the decay z/g —^ t'e + takes place, the period is too short to have signihcant energy transfer 
and therefore the energy loss due to z/g —> he + does not imply any bound on \gee\- In 
the next period (t > 1 sec) when Vj, < 0, neutrino decay is no longer kinematically allowed 
and instead the two processes Pg —z/g + J and z/g + z/g —J can take place. Since the 
density of z/g is much higher in the inner core, the process z/g + z/g —J implies a stronger 
bound. We have found that l^fggl < 10“^ (see Eq. (42)) if the emitted Majoron leaves the 
core immediately after production. We have found that for \g\ < 10“^ the effect of four- 
point processes (z/g -|- z/g ^ J -|- J and Pg -|- z/g —J -|- J) is negligible. We believe that 
previous treatments of the reaction [e.g., [15]) have not correctly subtracted the 3-particle 
subprocesses. 

We also have studied the bounds on coupling of Majoron to muon (tan) neutrino. In 
the basis in which g^^. = 0 (since z/^ and z/^ are equivalent for supernova processes, we can 
always rotate (z/.^, t'/i) to a new basis (z/j., z/^) for which gu'^v'^ = 0) we have found the following 
results. For \geij\'^l\g^j.^j\^ x < 37, the processes P/i —z/^ + J and z/^ + z/^ —^ J imply 

Id/i/il < 8 X 10“’^ {see Eq. (49)) while z/g -|- z/^ —J gives \ge^\ < 5 x 10“'^ {see Eq. (50)), 
providing the emitted Majoron leave the core without being trapped or undergoing decay. 
For \ge^i \‘^X {iil^/T"^ > 37, we have shown that the process z/g -|- z/^ —J eats up z/^ 
within a short period (leading to a negative chemical potential for z/^) such that the bound 
from z/g -|- z/^ —J no longer applies. However, in this case, the density of P^ increases 
~ becomes positive) and the bound on (Eq. (49)) still applies (actually 

it will be a conservative one). For \ge^\ > \g the P^-decay in the outer core 
(where g,u^/T ~ 1) imposes the strongest bound on \ge^i\ which is \ge^l\ < few x 10“®. These 
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upper bounds are schematically summarized in Figs. (1) and (2). Note that the bounds on 
\gTT\ and \geT\ are exactly the same as those on and \getj\i respectively. The bounds are 
parameterized in terms of supernova parameters (T, chemical potentials, V and the radius 
of the core) so it is possible to apply the predictions of any supernova model. 

All these upper bounds come from three-particle processes shown in diagrams (a,b) of Fig. 
(2). In these processes, all the involved particles are on-shell. Therefore the aforementioned 
bounds can be translated into bonnds on the corresponding element of the matrix h in the 
derivative form of interaction (see Eq. (2)), nsing the relation given in Eq. (6). 

We also stndied Majoron decay and the interactions that can trap Majorons. We have 
fonnd that the processes z/g + J —> ^ and z/^^^ + ^ niay have signihcant 

effect > Rcore)y Only if \gee\ > 6 X 10“®(g/10 MeV)^/^, \ge^l{T)\ > 2 x 10“®(g/10 MeV)^/^ 
and \g^(r)^J,{T) \ > 4 x 10“®, respectively. If the conplings of Majorons to nentrinos are larger 
than these limits, the Majorons cannot leave the core immediately. However, the processes 
involving the Majoron still affect the evolntion of snpernova, transferring energy from the 
inner core and distorting the density distribntion of the particles. If the conplings of Majoron 
are larger than some lower bonnds, the only Majoron particles that can leave the core and 
cool down it are those prodnced in (or diffnsed into) a shell close to the nentrino-sphere 
where the density decreases rapidly with increasing radins. In this region the density is too 
low to give rise to a signihcant Majoron hux (he., Cj Cy). We emphasize that to derive 
the lower bonnds, it is not snfhcient to consider the conpling constants collectively. For 
example, if \gefj\ > 5 x 10“®, the Majorons prodnced via z/e + z/g —> J can annihilate with 
another z/g into before escaping the core. 

To derive the lower bonnds, one must recalculate the density and temperatnre prohles of 
matter, nentrinos and Majoron particles which, in general, are different from those calculated 
so far without including Majoron processes. Here, we have discussed and evalnated only the 
fonr-point interactions which for large values of coupling constants may have signihcant 
ehect. 
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Figure 1: Bounds on g^e from supernova cooling (upper lines) and lab observations (lower 
lines). Line (a) shows the “upper” bound on gee from the z/gi/e —>■ J process in the supernova 
core (see Eq. (42)) while line (b) represents the “upper” bound from Ve ^ VeJ [14], Line (c) 
shows the “lower” bound which is derived without considering the effect of the four-point 
processes [14], Line (c) gives the “lower” bound according to [15], we have argued that this 
is an overestimation. Thus we expect the true “lower” bound (e) to be between lines (c) and 
(d). The range of parameters between the “upper” and “lower” bounds (the horizontally 
shaded area) is excluded by supernova considerations. Line (f) represents the upper bound 
from double beta decay [12] and the whole region to its right (the vertically shaded area) 
is excluded. Lines (g) and (h) represent the upper bounds derived from solar neutrinos [10] 
and Kaon decay [13], respectively. Note that the bound (f) from solar data applies to g 2 i 
rather than gee-, we have included this line to compare the orders of magnitudes of different 
bounds. We have not resolved whether the “lower” bound (e) lies above or below (f); in the 
latter case, there is a small allowed region between the two bounds. 
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Figure 2: Diagrams (a) and (b) are the dominant three-point processes and are possible only 
for V < 0. Diagrams (c) and (d) are the subdominant diagrams and can take place for any 
value of V. 
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Figure 3: The bounds on coupling constants for T = 10 MeV and = 200 MeV. The 
shaded area is excluded by energy loss considerations. The horizontal and vertical lines at 
5 X 10“® and 3 x 10“® represent the upper bounds obtained in Eq. (52) and Eq. (49), 
respectively. The dashed lines show the limits above which Majorons with energy ~ 10 MeV 
scatter before leaving the core. The dotted lines represent the same limits for Majorons with 
energy ~ 200 MeV (see Eqs. (62,63)). The dot-dashed line schematically represents the 
“lower” bound. We have not calculated the exact numerical value of the lower bound, but 
this is an estimate for = 0. Note that the energies of Majorons produced via —>■ J 

and ^ are of the order of 10 MeV; that is why the “lower” bound can be to the left 
of the vertical dotted line. The values inside the notch (above the horizontal line at 7 x 10“^) 
are excluded due to ^ J (see Eq. (50) and its discussion). 
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Figure 4: The bounds on coupling constants for T = 20 MeV and = 200 MeV. The 
shaded area is excluded by energy loss considerations. The horizontal and vertical lines at 
4.5 X 10“® and 8 x 10“^ represent the upper bounds obtained in Eq. (52) and Eq. (49), 
respectively. The dashed lines show the limits above which Majorons with energy ~ 10 MeV 
scatter before leaving the core. The dotted lines represent the same limits for Majorons with 
energy ~ 200 MeV (see Eqs. (62,63)). The dot-dashed line schematically represents the 
“lower” bound. We have not calculated the exact numerical value of the lower bound, but 
this is an estimate for gee = 0. Note that the energies of Majorons produced via —>• J 

and —>■ are of the order of 10 MeV; that is why the “lower” bound can be to the left 
of the vertical dotted line. 
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